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ABSTRACT 

The damping of primordial perturbations at small-scales gives rise to distortions of the cosmic 
microwave background (CMB). Here the dependence of the distortion on the different types 
of cosmological initial conditions is explored, covering adiabatic, baryon/cold dark matter 
isocurvature, neutrino density /velocity isocurvature modes and some mixtures. The radiation 
transfer functions for each mode are determined and then used to compute the dissipative 
heating rates and spectral distortion signatures, utilizing both analytic estimates and numer- 
ical results from the thermalization code CosmoTherm. Along the way, the early-time super 
horizon behavior for the resulting fluid modes is derived in conformal Newtonian gauge, and 
tight-coupling transfer function approximations are given. CMB spectral distortions caused 
by different perturbation modes can be estimated using simple A:-space window functions 
which are provided here. Neutrinos carry away some fraction of the primordial perturba- 
tion power, introducing an overall eflSciency factor that depends on the perturbation type. 
It is shown that future measurements of the CMB frequency spectrum have the potential 
to probe different perturbation modes at very small scales (corresponding to wavenumbers 
1 Mpc^' ^ k < few X 10"* Mpc"'). These constraints are complementary to those obtained at 
large scales and hence provide an exciting new window to early universe physics. 
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1 INTRODUCTION 



It is well-known that energy release in the early Universe leads 
to spectral dis tortions (SD s) of the CMB ( Zeldovich & SunvaevI 
, 1969; Stmyaev & Zeldovich 19701^; llllarionoy & Sunvaeyl 1 197a 



Danese & de Zottill 19821 ; iBurigana et al.lll99ll : IHu & Sil]dll993ah . 
At early times, a //-type distortion is created, while energy release 
at lower redshifts (z < 5 x 10'') results in a y-type distortion, sim- 
ilar to the Sunyaev-Zeldovich effect from galaxy clusters. Con- 
straints on the //- and y-parameter obt ained with COBE/FIRAS 
JMather et al.l 1 19941 ; iFixsen et aU 1 19961) limit possible deviations 
from a bla ckbody to n < 9 x 10"^ and j; < 1.5 x 10"' at 95% 
confidence l IFixsen et aljl996l) . At slightl y lower frequencies, there 
are similar limits to // from the ARCADE (Kogut et al. 2004.'200d; 
Seiffertetal]|201lh and TRIS jZannoni et al. 2008 ; Gervasi eT^ 
20081) experiments. 

Improvements in experimental design and technology may 
soon allow much more sensitive (a factor of ^ 10^ - 10* in 
/I beyond the FIR AS limit!) measurements of the CMB fre- 
quency spectrum, as suggested for the proposed experiment PIXIE 
JKogut et alj201 ID . This has spurred renewed t heoretical interest in 
the cosmological thermalization problem (e.g.. lChluba & SunvaevI 



* E-mail: jchluba@pha.jhu.edu 
t E-mail: dgrin@ias.edu 



I2OI2I ; IPaier & Zaldarriagal 1201 21; iKhatri & SunvaevI 1201 2bl) and 
the use of SDs to probe new physics. For instance, en- 
ergy injection from dark matter annihilation o r decay at red- 
shifts z < few X 10'^ JHu & Silkl Il993bl; iMcDonald et"d] 
I2OOII; IChlubal I2OIOI; I Chluba & SunvaevI l2012h. c osmic strings 
I Ostriker & Thompson Il987l; Tashiro et al.l l2012bl lal). primordial 
magnetic fields 1 Jedamzik et al. 200( j), but also more exot ic pos- 
sibilities jLochan et al.,.2012 ; .Bull & Kamionkowskill2013l ) could 
produce a detectable SD signal. 

It has long been known that Silk damping JSilM Il968l) 
of primordial small-scale pe rturbations also causes ener g y re- 
lease in the early universe dSunvaev & ZeldovichI Il970al ; iDalvl 
I1991I ; iBarrow & Co"ieslll99ll ; IHu et al.lll994al) . The SD signal de- 
pends on the amplitude and shape of the primordial fluctua- 
tion power spectrum, which in turn depend sensitively on the 
early-universe physics seeding these perturbations. The detec- 
tion and characterization of CMB SDs could thus offer a pow- 
erful n ew probe of inflationary models (Chluba & Sunyaev 20121; 
Khatri et al. 2012a; Paier & Zaldarriaga 2012; Dent et al. 20i: 



Ganc & Komatsi] |2012| ; IChluba et alj l2012d ; IPowelll |2012|) . The 



proper interpretation of measurements of the CMB frequency 
spectrum in the context of early universe physics requires an 
ac curate treatment of a coustic mode dissipation (as developed 
in IChluba et alj 1201 2bl ; iKhatri et alj l2012bf ; IPaier & Zaldarriag j 
|2013|) . which we follow here. 
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The amplitude o f the di stortion also depends on the type of 
perturbation mod es ( jB arrow & Coles 1991; Hu & Sugiyama 1994i 
IPent et al.ll2012r) . The combined system of fluid + Einstein equa- 
tions for cold dark matter (CDM), baryons, neutrinos, and photons|j 
has a variety of propagating (and growing) normal modes. 

Empirically, the most prominent is the adiabatic (AD) mode, 
for which the curvature perturbation f is non-zero, but all initial 
entropy fluctuations of other species relative to photons vanish: 
Siy - dUilrii - Snylriy = (5p;/[(l + Wi)pi\ - {3/4)Spy/py = 0. Here 
Wi denotes the equation of state of the ;"' species. The other well- 
behaved fundamental modes (on super-horizon scales) are known 
as isocun'ature or entropy fluctuations, and are characterized by 
initial values ^ = and 5 ,y i^ 0. They may be simply thought of as 
spatial fluctuations in the composition (or its time derivative) of the 
Universe (e.g., see ,Ma & Bertschinger 1995; Bucheretal. 2000, 
for additional details on perturbation equations and definitions). 
The most straightforward isocurvature modes are density isocur- 
vature fluctuations, which naturally sort into the CDM isocurva- 
ture (CI), the baryon isocurvature (BI), and the neutrino density 
isocurvature (NDI) mode. Additionally, there is a neutrino velocity 
isocurvature mode (NVI), with initial ^ = 6p = 0, but a relative 
velocity iv - i^ perturbation between photons and neutrinos. 

The standard lore is that isocurvature fluctuations behave 
differently from adiabatic perturbations, sourcing a much larger 
Sachs-Wolfe effect at large angular scales in the CMB, and acous- 
tic peak phases th at are out of phase with those produced by adia- 
batic fluctuations iIHu & Sugiyam3l 199511 19961 ; iKodama & Sasaki 
ll986l ; lEfstathiou & Bondll98dll987h . Indeed, this is true of BI and 
CI perturbations. NDI modes, on the other hand, behave more like 
adiabatic modes, since the initially perturbed species (neutrinos) is 
relativistic, and the isocurvature condition (6p = 0) requires an ini- 
tial energy density perturbation in the photons {Spy = -6py). Sim- 
ilarly, for the NVI a significant CMB dipole is present. The phase 
structure of the CMB anisotropy transfer functions for NDI and 
NVI perturbations thus bears a closer resemblance to the adiabatic 
transfer functi on than do the transfer functions for BI/CI modes (for 
examples, see lKawasaki et alj|2012r) . 

At large scales, corresponding to 10"* Mpc"' < k < I Mpc"', 
these important differences are readily probed by precise mea- 
surements of CM B anisotropies, an d are highly co nstrained by 
the B oomerang jjones et aU 1 20061) and WMAP JLarson et aU 
UOIlh expe riments, a s well as large-scale structure (L SS) mea- 
surements I Beltran et al. 20 05; Seliaketal. 2006; Zun ckel et alj 
I2OIII; lMuvaKasandaetalJl2012l) and other CMB experiments 
JMoodlev et alj |20()4| ; iMacTavish et al] I2OO6I ; iBean et alJ I2OO6I ; 
IPunklev et alj|2009f) . We thus know that at large scales the pri- 
mordial fluctuations are predominantly adiabatic, and there is no 
evidence of a significant isocurvature component to the primor- 
dial i nitial conditions. Recent results f rom the Planck collabora- 
tion jPlanck Collaboration et al.ll2013cr) stren gthen this view, al- 
though the hemispherical power asymmetry dEriksen et al.ll2004l ; 
IPlanck Collaboration et al.l2013dn could be interpreted as a hint for 
a modulated large-scale isocurvature mode (IDai et al.ll2013h . 

The implications of this fact are dramatic. In the simplest in- 
flationary models, fluctuations in all species are seeded by quan- 
tum fluctuations of a single scalar field; the consistency of ob- 
serv ations with adiabatic initial condition s supports these scenar- 
ios teardeenet all 19831; lOuth & PillT98^ ; lLvthlll984 iMaldacenal 



' In the following, for associated variables we use subscripts c, b, v and y, 
respectively. 



l2003h . Indeed, the Boomerang, WMAP and Planck data impose 
very strong constraints on topolo gical defect-domina t ed mod- 
els for cosmic structure formation ( Contaldi et al.lll999l ; lAlbrechtl 
I2OOOI ; lFraissell2005] ; IPlanck Collaboration et alj|2013d) . Therefore, 
a detection of a sub-dominant but non-zero isocurvature compo- 
nent would imply deviations from the simplest inflationary pic- 
ture. For example, dark matt er could be composed of axions 
l lSikiviell2003 ; ICadamurdl20 1 2|) . The axion field would be present 
and e nergetically subdominant during inflation, exciting the CI 
mode ilAxenidesetalJll983l ; lLindelll985l ; ISeckel & Tumeilll985l; 
[Turner & Wilczekll99lh . Isocurvature constraints thus limit the pa- 
rameter space available for axion dark matter, and SDs might help 
shed light on this. 

Alternatively, in the curvaton model, a sub-dominant scalar 
field (the curx'aton) picks up quantum fluctuations, comes to dom- 
inate the cosmic energy budget, and then seeds a correlated mix- 
ture of adiabatic and isocurvature fluctuations. The amplitudes of 
CI, BI, and NDI fluctuations are then set by the relative placement 
of the epochs of dark matter production, lepton number creation, 
or b aryon number creati on, relative to the time of curvaton de- 



cav JLvth & Wan ds 2002; Lvth et"ani2003l ; lG"ordon & Lewill2003l; 
[Gordon & Pritchardi2009.> . It would be useful to determine if SDs 
off'er any additional leverage on curvaton parameter space. 

The curvaton model may also excite compensated isocun'a- 
ture perturbations (CIPs), for which Sp = ^ = 0, and 6S - 
dn^/n^ - Sn^jjn^, + 0. Surprisingly, current CMB data analyses do 
not impose constraints to CIPs, but they coul d soon be detected 
using higher-orde r correlations of the CMB (IHolder et al.ll201Cl ; 
lOrin et al.ll2011allbl) . It would be interesting to see if SDs could 
be used to detect CIPs or other curvaton-induced modes on small 
scales, a possibility analyzed in this work. We find that SDs are 
only useful for this purpose at very futuristic sensitivity levels. 

Another interesting theoretical possibility is inhomo geneous 
baryogenesis, which predicts the existence of BI modes jPeeblej 
Il999j d; iKovama & s"od3 1 19991) . More broadly, multi-field in- 
flationary models e xcite isocurvature fluctuations (IGordonluOOll ; 
I Gordon et al.l uOOll) . Constraints to all these possibilities from 
CMB data are informative, and promise to be even more sensi- 
tive with the next Planck cosmology data release, including the full 
temperature and polarization information. 

CMB SDs could allow the characterization of primordial fluc- 
tuations on length scales far smaller than possible with CMB 
anisotropy measurements. Most recent work on CMB SDs explores 
their dependence on the power spectrum of primordial fluctuation s, 
restricting attention to the adiabatic m ode dChluba et al.ll2012bl lat 
|Powellll2()lj ; iKhatri & Sunvae-^j|2013l) . SDs could, however, also 
test for the presence of isocurvature initial conditions, providing 
a complementary probe to measurements of CMB anisotropies 
and cosmological large-scale structure on radically different length 
scales than those measurements. 

The imprint of isocurvature fluctuations on the CM B spec- 
trum was first studied in detail by IHu & Sugivamal ( Il994) . where 
BI models were tightly constrained. In BI models, stars would form 
early, leading to early reionization and a detectable Compton y- 
type distortion of the CMB frequency spectrum. More recently, the 
early energy rele ase from BI/CI fluctuations was investigated by 



i ' rele a? 
20l3), 



iDent et al.l i2012h . using a simplified treatment of both the trans- 
fer function that maps primordial fluctuations to moments of the 
radiation field, and of the heating rate itself. 

Recent work has shown the heating process is only 3/4 as 
efficient as previous estimates, and that closer to recombination, 
baryon loading and second order Doppler terms are important to the 
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heating rate calculation JChluba et al.ll2012bf) . Furthermore, free- 
streaming relativistic particles, like neutrinos, carry away some 
fraction of the perturbation power, introducing a dependence on the 
effective number of relativistic species. These details are carefully 
considered here but were omitted in previous work. 

Here we revisit the problem, using the pe rturbation and heat- 
ing modules of the SD c ode CosmoTheriv|j JChluba & SunvaevI 
I2OI2I : IChluba et alj|2012bl) to precisely calculate both the evolu- 
tion of fluid and radiation variables from isocurvature initial con- 
ditions, as well as the resulting effective plasma heating rate and 
CMB spectral distortion signal. We comput e the SD signal fro m BI 
and CI models, improving the estimates of lOent et al.l l l2012r) . and 
extend our reach to the SD signature of NDI and NVI modes, as 
well as the CIP ( baryon-CDM isocurvature) mode. We expand on 
iDent et al.l 1 120 12h to compute both yu and v-type SDs for adiabatic 
and isocurvature modes, exploring the dependence of the signal on 
the spectral index of the initial power spectrum of each mode. We 
compare these different possibilities with the sensitivity of the pro- 
posed PIXIE mission and existing limits to the chemical potential 
/J and Compton y-parameter from COBE/FIRAS. We also provide 
simple analytic expressions for the heating rate of different pertur- 
bation types, as well as A-space window functions which may be 
used to estimate the SD {ji and y) signal of arbitrary power spectra. 

We begin in Section |2] with a review of acoustic mode dis- 
sipation and the resulting plasma heating responsible for spectral 
distortions (SDs). In Section [3] we move on to discuss the differ- 
ent families of cosmic initial conditions, and the resulting acoustic 
mode amplitudes at small scales. We also consider mode mixtures 
in this section, but restrict ourselves to the simplest cases, provid- 
ing a simple recipe for correlated, uncorrelated and anti-correlated 
modes. In Section|4l we present the precise numerical heating rates 
produced by different acoustic modes. In Section|5]we discuss pos- 
sible constraints to the power spectra of different pure modes and 
some representative mixtures of different modes, as well as im- 
plications of SD experiments for specific early universe scenarios, 
such as the curvaton model. We conclude in Section|6] 



2 CMB SPECTRAL DISTORTIONS CAUSED BY THE 
DISSIPATION OF ACOUSTIC MODES 

In this section we briefly review how the dissipation of acous- 
tic modes c reates spectra l distort ions. For a more in-depth dis- 
cussion, see IChlubaetall l l2012bh and iKhatri et alj ( l2012d) . The 
problem boils down to computation of the effective heating rates 
for different perturbation modes. These are obtained both numer- 
ically, by solving the co smological perturbation equations (see 
iMa & BertschingenI 19951 for details) with CosmoTherm, and an- 
alytically in the tight-coupling approximation. With the heating 
rates in hand, we use simple analytic estimates to compute the 
resulting chemical potential, ^, and Compton y-parameter. Gen- 
erally, the detailed shape of the distortion is not just represented 
by a simple superpositio n of /i- and y-distortion, as shown in 
IChluba & Sunvaev (2012, e.g., see Fig. 15 and 19) and more re- 
cently by Khatri & Sunvaev (2012a). These details will, however, 
be addressed in some future work, since for estimates the approach 
presented here suffices. 



2.1 Estimates for the /j and y parameters caused by early 
energy release 

At high redshifts, z » z^ x 1.98 x 10^, the thermalization process 
is extremely efficient. As a result, any energy release just increases 
the specific entropy of the Universe, and thus raises the average 
temperature of the CMB without producing SDs. For lower red- 
shifts, z < z^, thermalization becomes less efficient, and energy 
release can produce SDs. For all z -^ Zfi, energy injection initially 
appears as a y-distortion to the CMB blackbody. If additionally, 
z 5S> Zp ,, a 5 X 10"*, Comptonization of the radiation field is still ef- 
ficient and the initial y-distortion is mostly converted into a chem- 
ical potential /j. On the other hand, if z ^ z^,,, Comptonization 
is inefficient, and the SD take the form of a non-zero Compton y- 
parameter, and essentially amounts to an early-universe analog to 
the Sunyaev-Zeldovich effect. Thus a yu-distortion is created by en- 
ergy releaseat5><20^ < z < 2x 10'' and ay-distortion at z < 5x lO'* 
(see lHu&Sii3ll993J. for more details). 

To estimate the values of the chemical potential, /j, and Comp- 
ton y-parameter, it is sufficient to compute the effective energy re- 
lease during the con'esponding epochs. Defining the distortion visi- 
bility function, SFhbiz) = exp (-[z/Z(i]^''), the weighted total energy 
release in the yu- and y-era is 
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We introduced the energy 
mordial acoustic modeij 



release caused by the dissipation of pri- 
a^*p7,^ d(a''(2ac)/dz, however, any pro- 
cess leading to energy release can be added here. Here the factors of 
the scale factor (normalized to unity today) a = (1 -1- z)"' cancel the 
main redshift dependence of the background radiation field, which 
is irrelevant for the creation of SDs. The factor Jbb(z) parametrizes 
the thermalization efficiency accounting for the effects of pho- 
ton production/d estruction by double C ompton sc attering (see 
Danese & de Zotti 1982; Burigana et al.|[r991: Hu & Silklll9933 : 



Chluba & SunvaevI 120 121 : IKhatri & SunvaevI I2012bl . for more de- 



tails). 

The simple expressions from ISunvaev & ZeldovichI ( Il970br) . 
fi X lAApy/py\f, and y x jAp-y/py\y, Eq. ([T]l can be used to es- 
timate the expected distortion at high frequencies. This imposes 
upper limits to any energy-releasing process in the early Universe. 
The relative ratio of the yu- and y-parameter in principle can be fur- 
ther used to distinguish different sources of early energy release via 
their redshift dependence, although this is possible only for spe- 
cific models of the thermal history. At low redshifts, aft er the re- 
combin ation epoch, many me chanisms [ e.g., re ionization JHu et alj 



Il994bf) ; supernova h eating ( Oh et al. 2003h: large-scale struc 



ture formation shocks ([Sunvaev & Zeldovichll972tlCen & Ostriken 
ll999l : lMiniati et alJuOOOl) : unresolved Sunyaev-Zeldovich clusters 
and the warm hot i ntergalactic mediu m ([Markevitch et al. 199l|; 
Ida Silva et al.ll2000l : IZhang et alj|2004h l give rise to large average 
y-type distortions. As a result, early universe processes which gen- 
erate a yu-type SD are more readily constrained with measurements 
of the CMB frequency spectrum, although useful bounds may still 
be derived from measurements of or limits on y-type distortions. 
One way to distinguish y-distortions from the pre-recombination 



^ CosmoTherm is available at jwww . Chluba . de/CosmoTherm| 
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epoch from those created at lat er stages m ight be the cosmo- 
logical recombination radiation (iChluba & Su nvaev 2006, 2009J; 
ISunvaev & Chluball2009l : IChluball2010h . but a more detailed dis- 
cussion is beyond the scope of this paper. To make further progress, 
we must specify the effective heating rate, Qao and its dependence 
on the cosmological initial conditions. 



2.2 The efTective heating rate from acoustic damping 

Small-scale perturbations of the photon temperatu re are completel y 
erased by shear viscosity and thermal conduction ( IWeinberg|2008h . 
These processes isotropize the photon-baryon fluid and lead to 
the mixing of blackbodies with slightly differe nt temperatures 
jZel'Dovich et alj|l972l : IChluba & Sunvaevll2004l) causing an in- 
crease of the local average photon temperature and a y-type SD. 
The spatially averaged SD source function, (Sx), directly depends 
on the amplitude and shape of the primordial perturbation power 
spectrum, Pi(k), as well as the detailed evolution of moments of 
the radiation field for given initial condi tions. It is determined by 
JChluba et al.ll2012bl : iKhatri et alj|2012bl) : 
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k^dk. 



2n^ 



■P,(k) 



(30, - vY 
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-^@2{@l + &l) + J](2e+i)@j 



(2) 



where @( and 0^ denote the photon temperature and polarization 
transfer functions and v the one for the baryon velocity. This source 
function can be computed accurately using the cosmological ther- 
malization code CosmoTherm. 

Given (<Sac), the required effective energy release rate caused 
by the damping of acoustic modes is determined by 



1 d|a*e,,| 4t(.Sac) 
Hd+zY 



afpy 



dz 



(3) 



where t = a-jN^c x AA x 10"'' (1 + zf sec"' denotes the rate of 
Thomson scattering and H x: 2.lx 10"^" (1 -I- z)^sec"' is the Hubble 
expansion ratqj. The factor of 4 arises because a y-distortion causes 
a change in the photon energy density by Apy ^ 4py. The factor t 
arises because the source function, (<Sac>, is defined with respect to 
the Thomson-scattering time-scale and l/[H(l + z)] is needed for 
the conversion to dz. 

Below we compute the effective heating rate for different ini- 
tial perturbation modes, with particular focus on AD, BI, CI, NDI, 
NVI modes (see Sect. |4] and Fig. [3](, and simple mode mixtures. 
The important differences are caused by the transfer functions and 
their relation to the initial power spectra, which can be understood 
using some simple analytic approximations deep into the radiation- 
dominated era, when photons and electrons are tightly coupled. 

2.3 Source term before recombination 

In this work we are particularly interested in energy release well 
before the recombination epoch (z > iC). At that time, the Uni- 
verse is still radiation-dominated and with small baryon loading 
R = 3pb/4py ~ 673(1 -I- z)"' < 7%. Also, photons and baryons 
are tightly coupled so that (30i - v) ^ 0. Furthermore, higher or- 
der temperature perturbations with f > 2 are negligible (Thomson 

■* The approximations for t and H are only valid at high redshifts, during 
the radiation-dominated era. 



scattering isotropizes the radiation field) and the dissipation physics 
is mainly determined by the quadr upole anisotropy. In thi s limiu, 
t02 ^ j^kSi and 0^ + 0^ ^ | 02 JHu & Sugivamjl 19961) . so that 



rk-dk 15 , I r dk , 16 , 

<^..>^J^/',(^)^0^«-,J^^^P,(^)^0?. 



(4) 



This result shows that an approximation for the source term can 
be obtained using analytic expressions for the CMB dipole trans- 
fer function. Inside the horizon, 0o and 0i ^ -d,,Qo/k have the 
generic form ( IHu & Sugivamalll996r) 



00^ 
01 : 
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(i+Ry* 

(i+Ryi* 



[Aik) cosikr,) + B(k) s,m{kr,)] e"*'/*D 
[A{k) sin(/tr,) - B{k) cos(ytrJ] e-'''"'o (5) 



where (cJcY = 1/[3(1 + R)] ~ 1/3 is the photon-baryon sound 
speed, rj = j c dt/a denotes conformal time and the damping scale, 
ko, is determined by 



d,k-^' 



2t 



R^ 
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(6) 



Furthermore, r^ is the sound horizon, and the WKB amplitudes A{k) 
and B(k) are determined by the initial condition (see Sect.|3}. 

For a given wavenumber k the source function oscillates 
rapidly, but for the net effect of many modes on the CMB spectrum 
we are only interested in time-averaged values. Squaring Eq. lO 
and averaging over many periods, interference terms vanish and by 
replacing sin^(fcrs) — > 1/2 and cos^ihs) —> 1/2 we obtain 



\k) + B^(k)\ e 



x\d,k-^ ^ §;_k'P.ik)C 

~ r- \'Pi(k)C^{k)d,e-^'''^''id\nk. 



-ik^ikl 



\k)e 



2k- /ki, 



(7) 



Here we defined C\k) = A^(k) + B-{k) and Ptik) = PPi(k)/(27T-). 
This expression can be used to estimate the effective heating rate 
caused by initial fluctuations, characterized by their power spectra 
and perturbation type (Sec.|4). Furthermore, since a separation of 
time- and scale-dependent terms is achieved, it is possible to de- 
fine A:-space window functions that can be pre-computed once the 
cosmology is fixed (Sect. 15. It . Differences between various per- 
turbation modes are then determined by the overall normalization, 
C^{k), specifying a mode-dependent heating efficiency. 



3 DEPENDENCE ON THE POWER SPECTRUM AND 
FLUID MODES AT SMALL SCALES 

As equations lO, l[3) and Q show, the energy release de- 
pends directly on the primordial perturbation power spectrum, 
Pjjk), which we character ize using the simple parameterization 
JKosowskv & Tumeilll995h , 



Pi(k) = 27T^k-'Ai(k/koY 



aHk/ko) 



(8) 



with amplitude Aj, spectral index ;i,, running «,,„„ = dnj dink, 
and pivot scale ko, which we set to ko = 0.002 Mpc"' . 

The definition of the power spectrum, Pj(k) (in particular, 
whether or not Pi{k) describes density fluctuations or fluctuations 

' Notice that in our definition 0^" = (2f + l)0f . 
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of a gauge-invariant variable like ^ or the entropy S ,-,,) may vary, 
depending on the perturbation type (see below). For adiabatic per- 
turbations, SD constraints for a variety of theoreti cally instruc- 
tive val ues of n, and n, p„n were recently discussed by IChluba et alj 
( l2012bl lal). in the context of inflationaiy theory. Here we focus on 
the dependence of the heating rate on different perturbation types, 
and provide simple expressions for the sub-horizon amplitudes, 
C^(k) = A^{k) + B^(k), defined by Eq. {T}, of small-scale modes. 
These allow comparing the heating rates of different perturbation 
types and also provide simple means for understanding the depen- 
dence on the spectral index (Sect.Q. 

To determine the coefficients, A{k) a nd B{k), of Eq. Q, w e 
can resort to analytic approximations (e.g.. lHu & Sugivamalll99q) . 
or simply solve the e volution equations for the fluid and metric vari- 
ables (see lMa & Ber tschinger 1995), for definitions) numerically for 
different initial conditions (see Appendix lAl for pure modes) to de- 
termine the main dependencies on scale and cosmology. The evo- 
lution of the potentials and their decay during horizon crossing af- 
fect the mode amplitudes in a non-trivial way, and so we use the 
latter approach to obtain a more accurate but simple description. 
Given approximations for A{k) and B{k), the heating rates for gen- 
eral mode mixtures at high redshifts can be constructed. Here we 
restrict ourselves to simple mode mi xtures, although more general 
cases could be of theoretical interest ( IMoodlev et al.ll2004l) . 



3.1 Adiabatic mode (AD) 

For adiabatic (isentropic) perturbations, we consider the power 
spectrum of curvature p erturbations, P,-{k)- The WKB mode am- 
plitude at small scales is l IHu & Sugivamalll99q) 



( lHu&Sugiv"am3ll996h 



4 
15 



(9) 



and 5 =: 0. Here Ry = Pylipy + Py) ~ 0.41 denotes the fractional 
contribution of massless neutrinos to the energy density of relativis- 
tic species. The term ARyjlS accounts for the correction caused by 
anisotropic stress in the neutrino fluid. It allows neutrinos to carry 
away some part of the perturbation power, without sourcing any 
CMB spectral distortion. Increasing the effective number of neutri- 
nos therefore decreases the net heating rate and SD. 

Although the initial temperature perturbation of the monopole 
is about 3 times smaller than A, decay of the potentials after hori- 
zon crossing boosts the mode amplitude to this larger value by 
gravitational forcing. Numerically we find a small admixture of 
the sine term, i.e., fi ^ 0, to the photon monopole transfer func- 
tion, Eq. l|5](, caused by the driving term. As shown in Fig. [T] the 
overall amplitude of the small-scale mode is well represented by 
|C|^(l+4i;„/15)-',but 



1- 0.338 iJ, 



B^ -7.16 X 10 --0.418i;„ 



(10) 



provide better approximations for the sine and cosine terms. At 
Ry = 0.41 this gives A ^ 0.86, B ^ -0.24, and C^ ^ 0.81. The 
dependence of A(A;) and B{k) on R,, is illustrated in Fig.|2] 



3.2 Baryon and CDM isocurvature mode (BI/CI) 

For BI and CI perturbations the situation is very different. In this 
case, Pi(k) is defined as the power spectrum of density pertur- 
bations, and the WKB mode amplitudes are roughly given by 



S=a 



V6 a, k^ 



4 a„ 



— 1 - — i?, 

15 



(11) 



and A ^ 0. Here k^^ ^ 9.46 x W-'ilJr Vl-i^yMpc"' ^ 9.56 x 
lO""* Mpc ' is the wavenumber of a mode crossing the horizon at 
matter-radiation equality, and i = {b, c) for baryons and cold dark 
matter, respectively. As this expression shows, for scales k « fceq 
the amplitude of photon temperature perturbations is suppressed. 
This occurs because these modes enter the horizon during radiation 
domination, when the gravitational sourcing of photon temperature 
perturbations by baryon/CDM density fluctuations is suppressed. In 
this case, the dissipative heating rate is larger at late times (during 
matter domination), when smaller k modes enter the horizon. 

From the WKB solution, we can also see that the heating rate 
decreases as the number of effective neutrino species increase^ 
Numerically we find that the total mode amplitude is represented 
sfightly better by replacing (1 - ARyll5) -> (1 + 2Ry/5r\ but the 
difference is only a few percent. For the WKB amplitudes we find 

A^--^-2i[4.79x 10-' + 0.195/?„l 



— -^ -J- [0.613 -0.235 Ry], 



(12) 



which gives A ^ -0.13(n,/n,„)(A:e,//t), B == -0.52(nja„,)(k^/k), 
and C^ ^ 0.2S(njnj-(k,Jkf for Ry = 0.41 (see Fig.|2j. 

Our calculations also show that, in contrast to the AD mode, 
the potentials do not decay as fast after entering the horizon. Conse- 
quently, the zero point of the monopole transfer fu nction is offset by 
- i//, in agreement with prev ious analysis (e.g., see lHu & Sugivamal 
ll995l : lKawasaki et alj|2012r) . For the heating rate this aspect is not 
important, since the local mono pole does not source any significant 
distortion JChlubaet al.l2012bl) . 



3.3 Neutrino density isocurvature mode (NDI) 

Unlike BI/CI modes, NDI modes begin with non-vanishing initial 
(super-horizon) potential perturbations if) and if/, which immediately 
source photon perturbations close to horizon crossing. Addition- 
ally, as a result of the isocurvature condition (6p = 0), the initial 
neutrino density perturbation Spy requires an equal but opposite 
photon energy density perturbation Spy = -Spy (see Appendix lAt. 
The non-zero photon energy density perturbation means that in 
principle both the sine and cosine parts of the monopole transfer 
function are excited. Thus, perturbations in the neutrino density im- 
mediately start oscillating with appreciable amplitude after horizon 
crossing, rendering the sine part subdominant. From our numeri- 
cal solutions we find that at small scales the total amplitude of the 
neutrino density isocurvature mode is well represented by 

2Ryl Ry 



\C\- 



(13) 



with Ry 



1 - Ry 



5{\+3Ryl5y 
0.59. Since small-scale NDI modes behave 



similar to AD modes, the SD caused by NDI modes (up to an over- 
all efficiency factor) is also expected to be comparable. The cosine 
and sine amplitudes of the monopole transfer function are 



A ^ -(Ry/Ry) [0.316 - 5.25 x W'^ R,] 
B ^ -(Ry/Ry) [0.267 - 0.297 Ry] , 



(14) 



For small changes, ARyfRy <K 1, the scaling is similar to that of the 
adiabatic modes. 
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Figure 1. Amplitude of cosine and sine parts in tlie WKB approximation of the sub-horizon monopole transfer functions, Eq. (|5]. The CI mode amplitudes are 
equal to the BI amplitudes times n^dm/^b- We numerically solved the fluid and potential equations for k = 10 Mpc" as representative example for small-scale 
modes varying i?,,. 



which gives A == -0.2, 5 ^ -0.1, and C^ == 0.052 (see Fig.©. 



3.4 Neutrino velocity isocurvature mode (NVI) 

NVI modes are excited by initial perturbations in the neutrino ve- 
locity, and in this case Pi{k) is a power spectrum of the fluid ex- 
pansion of the photon-neutrino relative velocity, 6y - 6y. Like NDI 
modes, NVI perturbations immediately start oscillating upon hori- 
zon entry, exciting sine and cosine terms comparable with relative 
amplitudes, depending on Ry. In the limit of ^^ <K 1 mainly the 
cosine part of 0q is excited, while for Ry ^ \'H is the sine part (see 
Fig.[TJ. We find that the total amplitude is well approximated by 



3.5 Mixture of difTerent perturbation modes 

Assuming that initial perturbations are created by a single field, the 
different pure modes can be excited in different proportions. With 
the expressions for A{k) and B{k) given above it is thus straightfor- 
ward to compute the heating efficiencies. The WKB amplitudes are 
simply given by coherent superposition 



A{k) ^ Y, ".(fcM,(fe) 
B{k)^Y,P,{k)B.{k), 



(17) 



\C\- 



3 yf3Ry/Ry 
5(l + 6Ry/5)' 



(15) 



The overall heating caused by NVI mode (like the NDI mode) is 
thus expected to generate SDs similar to those caused by the AD 
and NDI mode. For the WKB amplitudes we find 



A ^ -(Ry/Ry) [0.935 - 1.06 Ry] 

B =^ {RylRy) [0.349 + 0.369 i?„] , 



(16) 



yielding A ^ -0.35, B ^ 0.35, and C^ ^ 0.25. The heating effi- 
ciency for NVI is thus ^ 5 times larger than for the NDI. 



where the sum runs over AD, BI, CI, NDI, and NVI mode ampli- 
tudes. The mixing coefficient, a, and yS, are determined by the ini- 
tial conditions. This superposition respects the phase of the transfer 
functions and the effective small-scale mode amplitude and heating 
efficiency is directly determined by C-(k) ^ A^{k) + B^(k). For the 
special case of totally uncorrelated modes (created by statistically 
independent processes), one finds 

CLo,-, - Yj »?(^)^'(^) + Yj^j(k)Bj{k) ^ Yj C], (18) 

/ ; ; 

where destructive and constructive interference terms average out. 

One simple example for a mode mixture is the CIP, for which 

(without loss of generality) the relative ratio of BI and CI am- 
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Isocurvature modes 



AD + NDI. con-elated 
AD + NDI, anti-conelated 
AD + NVI, correlated 
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\) 



§ 0.4 




Figure 2. Effective mode amplitude, C' , for mixture of AD with either 
NDI or NVI modes. We chose Ry = 0.41 for which we have Aad = 0.86, 
Bad = -0.24, Andi = -0.2, Bndi = -0.1, Anvi = -0.35, and Bnvi = 0.35. 



plitude is Oqi = ySci = 1 and ^bi = j8bi = -^cdm/f^b. which 
gives A{k) =^ B(k) ^ 0, and thus very small net heating. We con- 
firmed this statement numerically, finding that the CIP heating rate 
is suppressed by at least two orders of magnitude relative to the CI 
mode, which is itself expected to give a very small distortion (see 
Sect. l5.2t . We therefore omit CIP modes below. 

As an additional simple example, we consider mixtures of AD 
with either NDI or NVI modes. We can parametrize the relative 
amplitudes as Oad = j8ad = a = const and a, = Pj = i — a, where 
j e (NDI, NVI) and a e [0, 1]. For fully correlated modes we find: 

Acorr - oAad + (1 - a)Aj 
Sco„==afiAD + (l-«)5j. (19) 

Similarly, for fully anti-correlated perturbations we have 

Aanti - a A AD - (1 - ffM; 

S.„ti - aSAD - (1 - a)Bj. (20) 

These expressions directly determine the effective mode ampli- 
tude, C^(fe), which we show in Fig. [2] for Ry = 0.41. For the cor- 
related AD and NDI mode the net amplitude has a minimum at 
a =^ 0.17 or an isocurvature-adiabatic ratio of a/ (I — a) ^ 4.8, 
but due to the phase difference of the AD and NDI transfer func- 
tions the amplitude does not vanish completely. Comparing with 
the heating rate of the uncorrelated mode this is C^^^ - (1/4)C„„;,q^.. 
On the other hand, the anti-coixelated mode heating rate is about 
^anti - l-^C'„„(.(,jj at a = 0.2, with constructive interference dom- 
inating. Similarly, for the AD plus NVI mode the net amplitude 
becomes small for a ^ 0.35 or o-/ (I - or) =: 1.9. Relative to the 
uncorrelated mode amplitude this is C^^,.,. ^ 0. 13Cj^„^o,.,.. For the AD 
with NVI anti-correlated mode case we again find constructive in- 
terference with C^^ij ^ 1.9C„„;.ojj at a =; 0.35. 

The examples above are just meant to illustrate the computa- 
tion of the effective mode amplitudes, C^, at small scales for differ- 
ent mode mixtures. In a similar way one can consider AD plus CI 
or BI modes or even modes with three perturbation types excited. 
A detailed analysis is beyond the scope of this paper. 



comparing with results obtained with the analytic approximations 
[Eqs. (O, (O and mode amplitudes of Sect. [3). To initialize confor- 
mal Newtonian (CN) gauge fluid and metric variables in CosmoTh- 
ERM for the different perturbation modes, we require super-horizon 
power-series solutions for each mode, as a function of conformal 
time ;;. Using a matrix normal-mode analysis followed by a gauge 
transform, we obtain these power-series solution^ in Appendix IaI 

In Figure |3] we show our numerical results for pure adiabatic, 
BI/CI, and NDI/NVI modes. We assumed a scale-invariant pertur- 
bation power spectrum with normalization A, = I to allow a com- 
parison of the different heating efficiencies. For the CI mode we 
also varied the spectral index, iij. A cosmology with fi„, = 0.26, 
Qj, = 0.044, Qk = 0, ^A = 0.74, h = 0.71, ^p = 0.24, 
A^eff = 3.04, Tq = 2.726 K and improved recombinati on history 
JChluba & Thomasll20I ll ; lAli-Haimoud & Hiratall20I ih was used 
in all cases. 

Since for scale-invariant curvature perturbations the effective 
heating rate caused by adiabatic perturbations scale s roughly as 
d(gac7Py)/dz oc (1 + z)'^ at early times (e.g., see iKhatri et alj 
2012aj), we multiplied all rates by (I + z). As Fig. [3] shows, for 
scale-invariant primordial fluctuations not only for AD but also for 
NDI/NVI modes this means (I -I- z)d(2ac/Py)/dz - const before 
recombination. The sub-horizon mode amplitudes for AD, NDI, 
and NVI modes only depend on Ry. As a result, C" - const [cf. 
Eqns. ^, ([nil, and ^]. Using Eqs. (Hjl and Q, the (1 + z)"' red- 
shift dependence of these heating rates can thus be easily derived. 
On the other hand, BI and CI modes show a steeply increasing 
heating rate towards lower redshift. This is because for n, == I the 
amplitude of the photon temperature perturbations is suppressed 
by - ^cq/fc [see Eq. UlU at early times, and C-(k) scales like 
= (feeq/^)^- Consequently, for iij = 3 (i.e., with strongly increased 
small-scale density perturbations) one again expects a nearly con- 
stant heating rate (1 -I- z)d(2ac/Py)/dz - const, as found in our 
calculation (cf. lower panel of Fig. [3}. 

The dependence of the pre-recombination heating rate on the 
type of perturbation can be captured by substituting the fits of Sec. [3] 
into Eq. (O and defining an effective spectral index n', yielding 



D' 



I 



(I+4i?,,/15)2 



" 8in,„A:o 

r\2 'J I ^^c "-eq 



D' 



D' 



0.81 



0.19 



8 ma,: 



ySv - 4.7 



: Hi ■ 



■■ rij - 2 



25(1 + 3Ry/5y 

21{RylRyf 



0.05 



0.23 



(AD) (21a) 



(BI) (21b) 



(CI) (21c) 



(NDI) (2 Id) 



(NVI) (21e) 



25(1 + 6Ryl5y- 

withySv = (1 - 4/?,,/ 15)^ - 0.79. At early times, we therefore have 
the effective heating rates 



1 d\a*Q.A 



a*Py 



dz 



^2D- 



i"' 



{k)d:e--'''"'Ddlnk 



(22) 



where 'P*(yt) = A, (k/ko)"'-'^^-"--'-' '"**/*"'. We confirmed numerically 
that at redshifts z ^ 10"' these approximations work pretty well, 
giving ^ 10% - 15% precision for the effective heating rate. Here 



4 HEATING RATES FOR PURE MODES 

We now compute the effective heating rates for different perturba- 
tion modes using CosmoTherm, evaluating Eqs. ^ and ([3]l, and 



^ To our knowledge, this is the first published summai7 of early-time CN 
gauge super-horizon solutions for this system which includes all fluid vari- 
ables, in addition to the metric perturbations. 



© 0000 RAS, MNRAS 000, 000-000 



8 Chluba and Grin 



\i-y transition 



>• - era 



recombination 



10 F 



^ 10 



O) 



10 E 



a 10' 



10" 



10 = 



M .5 

S 10 



ta 



10 



10 







1 III 


; \3:y. 

r n.= l 

-_ I 

adiabatic 

— baryon isocurvature 


A -^ 

xf-ine ' 

— ■ • 

/ 
/ 

y 

y' 

y 1 


neutrino density 
— isocurvatme 

: neutrino velocity 

: isocurvature 

■ • 

1 1 1 1 1 1 ^ 


y 


.-' / 

/ / 

y^ 
y '' 

y X 

1 >i 1 



10 



10"^ 



a = [1+z] 



10" 







\i - era |i->' transition >• - era r 


ecom 


bination 


1 


- 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 IIIII 


'\ «i,„ = 3^ 


II 10 




I X/c~5.6 






O 








— ^ 


\ 


^ 10 


— 








\ — 


r-J 












"Td 




































e 10-' 


^ 






-■^' 


-^^-^^ 














9 10-^ 


: 


adiabatic, n = 1 

CDM isocurvature, 


.-■^'^^ 


/ 


^'\ 


•a ^" 




n. =1,2,3 


^.^" T^-^^o 


/ 


'^■^ E 




'- 


ISO 


^ '^ 


y 


' ■* ■ ^ '- 


- 




^•■^ •' 


' ^ 


^ 10- 


~ 


^."■^ 


i'' 




~ 


OJ 




^■^ 


■^^\ 




: 








y -^ 






2 .4 
M 10 


^ 


,.,.'-'~'(Uz)-^" 


y'^-^ 




^ 


a 


















y 


















ra 












OJ .> 












■c 10 


T^" 




y 




-= 


OJ 












> 






,/~[\^7f 




: 


O -6 












.iJ 10 






y 




-^ 


^K 












W 


- 


y 






- 


.^-7 


1 1 1 


1 1 -^ i 1 1 II 


1 1 1 i 1 1 III 



10 



10' 



a = [1+z] 



10" 



Figure 3. Effective heating rate, dlQjc/Pyl/dz, for different pure perturbation modes. We multiplied by (1 +z) and set the overall amplitude of the perturbation 
power spectrum to unity, i.e. A; = 1. We also used «;.,„„ = in all cases. For the upper panel we used spectral index n, = 1 while in the lower we 
varied it for the CDM isocurvature modes as labeled. The annotated factors are roughly /^b - (^c/^b)", /c - (3/8)(nc/n„,)^(^eq/<:o)"[l - (4R,,/15)^]-, 
/v - 4{Ry/Ry)^{l + 4i?,/15)-/[25(l + 3Ry/5)^] and /j =^ 27(1 + 3i?J5)2/[4(l + 6i?,,/5)-]. These illustrate the relative heating efficiencies for different 
perturbation modes with the same overall amplitude. The results were obtained by direct integration of the perturbation equations using CosmoTherm. 



e = 2D^ =^ const defines a mode-dependent heating efficiency. This 
implies that the early SDs produced by the different modes con- 
sidered here are all degenerate with an overall normalization when 
comparing AD, NDI, and NVI for n* = «, on one hand, with BI and 
CI modes for n* = «; - 2 on the other. The differences derive from 
how much of the initial perturbations in the different fluid variables 
at small scales actually appear as perturbations in the photon field. 
Comparing the heating efficiencies, Eq. l l21b . shows that AD 
modes dissipate their energy roughly 16 times more efficiently than 
NDI fluctuations. Similarly, NVI modes have ^ 4.7 times higher 
heating efficiency than NDI modes. Furthermore, BI modes source 
early SDs at about (Cl^/Clf,)- ^ 24 lower efficiency than CI modes, 
while in comparisons to AD modes CI fluctuations for n* = «, - 2 



cause ==5.6 times larger heating. All these statements are confirmed 
by our numerical results (cf. Fig.|3). 

Closer to the recombination epoch baryon loading no longer 
is negligible and we see a suppression of the heating rate relative to 
the high redshift scaling (cf. Fig.|3}- After the recombination epoch 
(z ^ 1000) the effective heating rates drop significantly as photons 
begin free streaming. At this late stag e the second order Do ppler ef- 
fect starts contributing significantly ( IChluba et al.ll20I2bh . For the 
baryon and CDM isocurvature modes the post-recombination heat- 
ing rate is relatively larger than for the adiabatic case, emphasizing 
the importance of late bulk motions for these perturbation modes. 
The neutrino isocurvature modes also result in relatively less heat- 
ing at the late stages. This implies that the relative ratio of /j and 
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y parameter is slightly mode-dependent, however, from the obser- 
vational point of view the differences are too small and degenerate 
with the shape of the power spectrum itself to allow distinguishing 
different scenarios in a model-independent way. 



4.1 Approximations for power law perturbation spectra 

One can gain additional insight by considering primordial fluctua- 
tions with power law perturbation power spectra («nin,i = 0). From 
Eq. l l22b it is straightforward to sh ow that the pre-recom bination 
heating rate is roughly given by (cf . IChluba et al.ll2012bl for adia- 
batic modes): 



1 difl^e,. 



a'^py 



3AiD^ /! + "* 



(1 + zf 



2klAo 



dz \+z \ 2 

where D- and ;i* > - 1 are mode-dependent [see Eq. ( 12 IH and 
(16/15)c 



(23) 



Ad«= 



18/fo^"^A'e,oo-T 



a; 5.92 X lO'^Mpc^ 



Here Hq is the Hubble parameter, fir is the density of relativistic 
species, and N^ = A'e.oCl + z)' is the number density of electrons 
(bound and free). Equation ( I23t explicitly shows that for AD, BI, 
CI, NDI, and NVI modes and «* = 1 the heating rate indeed scales 
like d((2ac/Py)/dz oc (1 -i- z)"'. In terms of the power spectrum, 
Pi{k), this mean scale-invariant perturbations for the AD, NDI, and 
NVI modes, while for BI and CI mode a very blue primordial power 
spectrum with n, = 3 is necessary. Also, for BI and CI fluctuations 
with spectral index n, == 1 («* ^ -I), from Eq. ([23} we find a 
redshift scaling diQac/Py)/ dz oc (1 -i- zy"^, while for ;?; = 2 («* = 0) 
one has diQ-^c/py)/ dz oc (1 + z)"^'^. Our numerical calculations 
confirm this dependence (cf. Fig. |3}, these scalings are, however, 
only valid at early time, before the recombination era (z > 10*). 
Especially when ;?, ^ 1, large scale modes contribute strongly to 
the total heating integral. The transfer functions of these modes are 
not well represented by the simple approximations given above, so 
that the integral formally diverges, unless a cutoff is introduced at 
small k. In this case one has to resort to full numerical integration 
of the perturbation equations. Still for iij > 1, Eq. ( 123b provides a 
fairly accurate estimate for the effective heating rate. 

Equation ( I23t also shows that the main dependence of the 
heating rate (for «;.,„„ = 0) for the adiabatic and neutrino isocur- 
vature modes is due to the overall amplitude of the perturbations 
power spectrum, Aj, the spectral index, n,, and the value of Ry, 
which changes the heating efficiency, e ^ 2D-. The dependence on 
the dissipation scale (related to A^ == 4.1 x 10"* Mpc"') is much 
weaker unless the small-scale spectral index, «,, differs strongly 
from unity. For the BI and CI modes Clt, and Q^ as well as feeq be- 
come important. Cosmology-dependence introduced by the dissi- 
pation scale is again less important for n, ^ 3. 



5 SPECTRAL DISTORTION CONSTRAINTS ON 
EARLY-UNIVERSE COSMOLOGY 

At large scales, k < IMpc"', constraints derived from CMB 
anisotropics and large-scale structure measurement are pretty 
tight already, suggesting nearly scale-invariant adiabatic pertur- 
bation with amplitude of the primordial cu r vature power spec- 
trum Ar ^ 2.4 x 10" (iDunklev et ai]|201ll: [KeTsler et alj[2oTII: 



IPlanck Collaboration et al.l l20I3ah . On th e other hand, at small 
scales constraints are much weaker (see iBringmann et al.l 1201 21 



for some discussion) and SDs provide a complementary (if not the 
only), strong probe for the primordial power spectrum. 

The discussion of Sect. [4] already indicates that the heating 
caused by very different perturbations modes and their mixtures, 
can lead to similar distortions. For example, given Pi(k) the time- 
dependence of the heating rates for AD, NDI, and NVI modes im- 
plies that from the practical point of view the SD signal should be 
indistinguishable, up to an overall efficiency factor that is degener- 
ate with the power spectrum amplitude. Similarly, the SD arising 
from dissipation of BI and CI modes should be practically indistin- 
guishable (small differences might arise in the post-recombination 
epoch, where non-linear effects will become important). Compar- 
ing the SD from AD, NDI, and NVI modes with those from BI and 
CI modes for given primordial power spectrum the former produce 
a much larger p distortion due to the ^ k^-^/k suppression of the 
BI/CI mode amplitude. This is, however, degenerate with the over- 
all spectral index of the perturbations. Still, on a model-by-model 
basis interesting constraints on the small-scale power spectrum can 
be derived for each case, as we illustrate here. 



5.1 Estimates for the p and v-parameters and the definition 
of A;-space window functions for pure modes 

We already showed that the pre-recombination heating rates from 
different perturbations can all be represented by one single ex- 
pression, Eq. ( I22t . with heating efficiency, e ^ 2D*, according to 
Eq. ( I21t . For power-law power spectrum one obtains the compact 
expression, Eq. ( |23b . Here we generalize to arbitrary shapes of the 
primordial power spectrum and estimate the SD by computing the 
redshift integrals, Eq. ([!). For the y-parameter the integral can be 
performed analytically, while for the /j-parameter the spectral dis- 
tortion visibility function, Jbb(z), requires numerical integration. 
This reduces the problem to a 1 -dimensional integral over A:-space 
window functions for the effective p and y-parameters: 



r°° k-dk 



{k) lyf (fc) 



'-^P,{k)W](k), 
where the A:-space window functions are 



iyf(A-)«2.8C'(A-) r J^^(z)d-j 



(24) 



(25a) 



~ 2.8 C^ik) 
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exp 



i 1360 J 



i + rxr + j , 

L 260 J 340 /' 



■ exp 



W]{k) 



, cHk) ^2^1^,,,^, ,_, c\k) 



exp 



(25b) 



with k = fc/[l Mpc" ] and cut-off scale, /;„]„ ^ 1 Mpc"', both be- 
cause modes at fc < I Mpc"' are already tightly constrained by 
CMB measurements at large scales, and because the analytic ap- 
proximations for the photon transfer functions introduced above 
become inaccurate. The approximations for Wik) are for the con- 
co rdance cosmo l ogy. Th ese expressions are similar to those given 
by IChluba et al.l ( l20I2ar) . but for W^ik.) we matched the numeri- 
cal results well not only in the small and large k limit but also at 
intermediate scales. In Eq. ( |25t , the fc-dependent factor C-(k) de- 
pends on the perturbation type. For AD and NDI/NVI modes one 
has C^(k) = D^ while for BI/CI modes C^{k) = D-(k) (ko/kf with 
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the values of D^ given by Eq. | |2U . The values of// and j-parameter 
for different modes can thus be estimated knowing D^ and rescal- 
ing the primordial power spectrum by appropriate powers of ko/k. 
Mode mixtures can be treated in a similar way (see Sect. 15. 3) . 



5.2 Constraints on different pure perturbation modes 

In this section we highlight constraints derived from SD measure- 
ments of COBE/FIRAS and future PIXlE-type experiments, focus- 
ing the discussion on pure perturbation modes. The derived limits 
should be interpreted as consen'citive upper bounds, since not only 
can several types of perturbation modes be present at small scales, 
but also other sources of early energy release (e.g., decaying or an- 
nihilating relicts, superconducting cosmic strings) could increase 
the CMB distortion. This would generally tighten the constraint on 
each source of energy injection. 



5.2.1 Pure AD, NDl, and NVI modes 

For AD modes a detailed discussion of SD power spectrum con- 
straints derived from COBE/FIRAS and a PIXIE-type experiment 
can be found in lChluba et alj l l2012ah . Since the pre-recombination 
heating rate for NDI and NVI modes only differs by an overall 
efficiency factor from the one of AD modes, their analysis di- 
rectly carries over For the NDI mode the relative heating effi- 
ciency is /„ ^ 4{RJR^fa + 4ie,,/15)V[25(I + 3RJ5f] ^ 0.061, 
which c an be captured by rep lacing A^ — » A^-//r - 16A^ in the 
work of IChluba et al.l J2012ah . Similarly, for the NVI mode we 
have / ^ 21{RJRy)-(\ + ARyjlSf l[25{\ + eRy/Sf] ^ 0.29 and 
A( — > 3.5 A(. Both for NDI and NVI modes the constraints are thus 
weaker than for the AD mode. 

To give some examples, for adiabatic modes PIXIE is able to 
rule out a scale-invariant curvature power spectrum with A^ > 4.3 x 
10"' at wavenumber k > 50Mpc~' with 5o '-confidence if no yu-type 
distortion is detected JChluba et alj2012bl lah. This therefore means 
that PIXIE would also be able to rule out scale-invariant neutrino 
density perturbations with overall amplitude A, > 7.0 x 10"* at 
wavenumber k > 50Mpc"' with 5(T-confidence. Even the yu-limit 
from COBE/FIRAS already implies A, < 1.3x 10"'' at 2cr-level. For 
the NVI mode these amplitude constraints are ^ 4.7 times tighter. 



5.2.2 Pure BI and CI modes 

Also for BI and CI modes the analysis of lChluba et alj ( l2012af) can 
be directly applied, however, not only are overall factors of /b - 
0.23 and /c ^ 5.6 needed, respectively, for the conversion of differ- 
ent constraints on adiabatic modes, but also one must use the effec- 
tive spectral index ;i* = n, - 2. For instance, for n' = 1 PIXIE would 
be sensitive to Ab > 1.9 x 10"* for the baryon, and A^ > 7.3 x 10"'° 
for the CDM isocurvature modes at wavenumber k > 50Mpc"' 
with 5iT-confidence. The limits from COBE/FIRAS are about 1800 
times weaker at 2cr-level. This value of «* means a very blue small- 
scale perturbation spectrum with spectral index n, = 3. This blue 
spectrum can be realized in certain axion isocurvature models, for 
example, in which the Peccei-Qui rm symmetry breaking sca le is 
itself dy namical du r ing in flation JKasuva & Kawasakil |2009| ), as 
noted in iDent et al.l l l20I2h . For «, = 1 the constraints are much 
weaker, as shown below. 



5.2.3 Simple expressions for power-law perturbation spectra 

For pure power-law primordial spectra, f, = A, (k/ko)"'^', we can 
simplify the computation of limits from p and y significantly us- 
ing the A;-space window functions, Eq. | |24K Defining the heating 
integrals 



U 



"^^I 



(yt//to)""' Wf(fc)dlnfc 



/,(«)= I {k/ko)"'^Wj{k)dlnk, 



(26) 



with Wf-{k) 



wf-{k)/D' and /t„ 



^ and kyniR = 1 Mpc"', limits on the over- 
all amplitude of the power spectrum at small-scales (k > I Mpc"') 
derived from p and y distortions can be expressed as 



A,< 



y"lii 



DH,{n',y 



A, £ 



y\u 



D%in;y 



(27) 



The y-limit probes power at scales I Mpc"' < k < 50 Mpc"', while 
the /j-limit is most sensitive to scales 50Mpc"' < k < W* Mpc"'. 
Late energy release (z £ lO"*), during and past the recombination 
epoch, is not included here. The corresponding y-distortion can, 
however, be computed using CosmoTherm in that case. 

In Fig. |4]we illustrate the dependence of /^(«*) and /,(«*) on 
the spectral index. These functions have a very step dependence on 
n*, which in the range -2 < «* < 5 can be approximated by 

ln/^(n) « 2.73 [l + 4.42^ + 0.444^^ - 9.21 x 10" Y' 

-0.0168f* - 5.38 X 10"^f^ + 4.92 x lO"**^] 

Inly(n) a 0.504 [l + 15.66^ + 0.845^^ + 0.0253f'' - 0.0189^] , 

with ^ = n - I. These expressions provide a 5% - 10% fit to the 
numerical results for Eq. (I26t . For « = I we find /^^(I) ~ 14.4 and 
/y(l) a 1.59 numerically. 

We can directly check the precision of these approximations 
for the heating integrals using CosmoTherm. The results are also 
shown in Fig. [4] For the /j-integral the approximations work very 
well. The approximations for the y-integral represent the full nu- 
merical result well if dissipation of modes at k < \ Mpc"' is ne- 
glected (consistent with the approximations made above). At larger 
scales, the approximations for the transfer functions are not valid, 
since the tight-coupling limit breaks down. Also, baryon loading 
and the second order Doppler effect become important so that we 
can only compute the effect accurately using CosmoTherm. Our 
result shows that a significant amount of energy is dissipated at 
late times if /i, < 1. These results are especially important for 
the baryon and CDM isocurvature modes, which for scale-invariant 
primordial density power spectrum would have n* = - 1 . This ren- 
ders a constraint derived from the y-parameter about ^ 1500 times 
tighter than the limit obtained by only accounting for modes with 
k > I Mpc"' and assuming that the power-law spectrum goes all 
the way from CMB scales to small-scales. These limits are still not 
competitive with those derived at CMB scales. For instance, from 
COBE/FIRAS we have |y| < 1.5 x 10"^ For n^ = -I this means 
Ab < 5 X 10"^ for the baryon and A^ < 2 x 10"^ at 2cr confidence. 
With PIXIE this could be improved by a factor of == 1800, but this 
is still far away from current constraints. 

For power spectra with n* > 1, constraints from p are gen- 
erally expected to be tighter than those obtained with y. Also, 
since there are many sources of large y-distortions at low red- 
shift, a clean detection of the primordial signal will be very chal- 
lenging. In all these cases one should consider the possibilities 
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Figure 4. Dependence of the heating integrals /,,(«) and /,(«) on the spectral index, n. We assumed pivot scale, ko = 0.002 Mpc" , values for k'^ ^ ko can be 
obtained by rescaling with (ig/to)'"". For comparison we also give the results for the heating integrals obtained with CosmoTherm for the adiabatic modes. 
We confirmed that the integrals for all the different modes considered here agree with each other to high precision. 



to include spacial-spectral and polarization information to sepa- 
rate the differen t components. For example, the pa tchiness of the 
y-distortion (see lPlanck Collaboration et al.ll2013bl for the first all 
sky y-ma p) from reionizatio n depends on the characteristic bub- 
ble size JZhang et al.l 12004) . which might allow primordial and 
reionization-induced y-distortions to be distinguished. A detailed 
discussion of this possibility is beyond the scope of this paper. 

Finally, in Fig. [5] we illustrate Icr constraints derived from 
COBE/FIRAS and PIXIE in the A,-w, plane. These were computed 
using Eq. ( I27t and JAU . Putting experimental obstacles aside, for 
AD, NDI, and NVI at n, > 0.9 the limits from /j distortions are 
tighter than those derived from a measurement of the y-parameter. 
For BI and CI modes this transition occurs at «, ^ 2.9. Further- 
more, due to the steep dependence on «;, constraints on the small- 
scale amplitude of AD, NDI, and NVI modes become very tight at 
«; > 1, showing the impressive potential of testing early universe 
models that produce large excess small-scale power. 



5.3 Constraints on mixed modes 

For the discussion in the previous section we assumed that only one 
type of perturbation was present. Constraints from CMB distortions 
on the amplitude of different modes of course can only limit the 
total energy release by the mode mixture. The heating efficiencies 
can be obtained as explained in Sect. 13.51 If correlations between 
modes are neglected, the final heating rates can be described using 
Eq. ( I22t and adding up the contributions of different modes. 

Mode mixtures can in principle give rise to interesting behav- 
ior. For example, a CI/BI isocurvature mode with very blue spectral 
index could be mixed with a quasi-scale invariant AD mode. This 
means that at small scales the isocurvature mode could dominate, 
while at large scales, which are well constrained by CMB experi- 
ments, the AD mode is most important. 

A similar story could be told for the NDI or NVI modes mixed 



with AD modes, and the constraints for these cases can be directly 
deduced from Fig. \5\ If on the other hand the isocurvature mode 
is completely subdominant, the SD limit will be very weak and 
mainly constrain the AD mode amplitude at small scales. For mix- 
tures with comparable contributions of AD and another mode, in- 
terference terms (destructive and constructive) make things more 
interesting. For example, in the curvaton scenario, isocurvature 
modes are all correlated with the adiabatic mode, and the precise 
coiTelation coefficients are determined by the curvaton energy den- 
sity fraction tq at curvaton decay, the lepton asymmetry parameter 
(which sets the number of effective relativistic degrees of freedom 
A'eff), and whether dark matter density, lepton number, and baryon 
number are produced before, by, or after curvaton decay. 

To determine whether or not curvaton-generated SDs could 
ever be detected, we sweep through the model space of 27 (3'') 
permutations, allowing the neutrino asymmetry parameter to span 
the entire (rather permissive) allowed experimental range . We gen- 
eralize Eqs. ([T9]l to allow all possible correlated mixtures of adi- 
abatic and isocurvature modes, and use the WKB coefficients of 
Sec. [3] We self-consistently compute th ese in the presence of a lep- 
ton asymmetry, a pplying expressions in lLvth & W ands (2002|) and 
iLvth et al.l 1 120031) . To properly compute the curvaton-generated y 
distortion, we use the numerically determined /,(«*)• 

The curvaton model would also seed local-type non- 
Gaussianity with amplitude fl^"\ We thus restrict consideration to 
the rang e rp > 0.15, the range still allowe d by Planck constraints 
to fl^l"' jPlanck Collaboration et al.ll2013eh . Additional limits to rjy 
could be obtained more directly from Planck constraints to isocur- 
vature modes, but these constraints have not been established in 
the p resence of all 4 isocurvature mod es with general correla- 
tions jPlanck Collaboration et al.ll2013cn . Even though the curva- 
ton model can excited the NDI mode, the overall change in /i and 
y over the null (adiabatic) hypothesis is of order 10%, and thus un- 
detectable at the sensitivity level possible with PIXIE. Future ad- 
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Figure 5. Limits on the amplitude of tlie power spectrum at ^ > 1 Mpc" for different pure perturbations modes and spectral indices. The heavy lines show 
constraints for a PIXIE-type experiment with lo" detection limits v = 2 X lO"' and fi = 10"**. Light hnes are present limits from COBE/FIRAS. Mode 
amplitudes above the corresponding lines are/will be ruled out. Assuming one overall power law perturbation spectrum at small scales, the limits derived from 
/i and y are not independent, and their ratio can in principle be used to distinguish AD, NDI and NVI on the one side from BI and CI, on the other (see text 
for discussion). Interpreting the limits independently, v distortions constraint power at 1 Mpc"' < k < SOMpc"', while the limit from /i probes power at 
50Mpc"' <k< 11000 Mpc"'. 



vances may change this dim state of affairs. An in-depth discussion 
of other correlated models is beyond the scope of this paper. 



6 CONCLUSIONS 

In the future, spectral distortions of the CMB might provide a pow- 
erful new probe of early-universe physics. Here we studied dis- 
tortions produced by the dissipation of small-scale perturbations, 
exploring the dependence of the signal on the different types of 



cosmological initial conditions. As one main result we obtained a 
unified formalism for the specific heating rates of the modes, allow- 
ing to describe the effect of pure modes but also mode mixtures in 
a quasi-analytic manner (see Sect.[3]and|4j. Our expressions can be 
used for precise computations of the SD signal using CosmoTherm 
to make more detailed forecasts, although here we restrict our at- 
tention to estimates of the associated /j and y-parameter, providing 
a simple way to constrain different early-universe models. 

We find that for scale-invariant initial conditions of compa- 
rable perturbation amplitude, the heating rates from pure BI, CI, 
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and CIP fluctuations are e xtremely subdomi nant to those from AD 
modes. In agreement with lDent et alj (20121) , we show that the SD 
signal of CI and BI modes falls below PIXIE's sensitivity, unless 
primordial perturbations have a very blue spectral index n, > 3. 
The BI SD signature is suppressed by ^ (flb/fic)" compared to the 
CI signal, a factor that is degenerate with the overall amplitude of 
the power spectrum. The CIP SD signal is strongly suppressed in 
addition by cancellation, and thus is unlikely to ever be detected. 
This all can be understood from the fact that neither baryons nor 
CDM can dramatically drive the evolution of the other fluid com- 
ponents until after matter-radiation equality. 

The NDI and NVI modes, on the other hand, yield heating 
rates comparable to the AD mode at all times, as neutrinos are rel- 
ativistic and similarly important to photons in driving the dynam- 
ics of the full coupled fluid system during radiation domination. 
We determined mode-dependent heating efficiencies [see Eq. J21b 1 
which weakly depend on the effective number of neutrino species, 
because additional non-interacting relativistic degrees of freedom 
carry away part of the initial perturbation power, never sourcing any 
perturbations in the photon fluid. Thus SDs are in principle sensi- 
tive to the presence of dark radiation or sterile neutrinos, however, 
this only causes an overall normalization factor that is degenerate 
with the power spectrum amplitude. SDs could in principle be used 
to probe parameters of the curvaton scenario for primordial fluctu- 
ations, although this would require rather futuristic 10%-level pre- 
cision in measurements of /j and y, which themselves are expected 
at the level A/// ^ 10"' - 10"* and thus are challenging to detect. 

While rather stringent limits can be derived for specific mod- 
els of the small-scale power spectrum (Fig. |5), SDs cannot tell the 
signature of different perturbation modes apart. For example, AD, 
NDI, and NVI modes should all cause a very similar SD signal up 
to an overall efficiency factor, assuming that the power spectrum 
has the same shape. This is because the heating rate at different 
redshifts scales in practically the same way (Fig.|3). Still, limits de- 
rived from SDs can be used as powerful tool to rule out different 
early universe models, and more detailed forecasts will be neces- 
sary to demonstrate the full potential of this new window. 
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APPENDIX A: INITIAL CONDITIONS IN THE 
CONFORMAL NEWTONIAN GAUGE 

CosmoTherm uses the confomial Newtonian (CN) gauge to follow 
the evolution of metri c and fl uid perturbations. U sing t he conventions 
of iMa & BertschingeJ j 19951) . iBucher et all J2000I) . and IShaw & Lewii 
)2O10l) , we can determine the required initial conditions for the different 



perturbation modes in this gauge, working deep in the radiation domi- 
nated epoch, and in the super-horizon regime. For convenience we define 
R,. = njnm,Rb = nb/nm,«r = \-Ry,ay = \ +4R,./15, &„ = 1 +2Ry/5, 
Cy = \ ~ ARyl\5, <j> = nn,//o/[4 Vni^c], and t = LdTj. According to these 
conventions, the scale factor evolves as (jJ aij) = r + t at times when only 
matter and radiation are energetically relevant. In these units, the scale fac- 
tor at matter-radiation equality is Ci> Oeq = 1/4. The overall normalization of 
the scale factor is iiTelevant to the final power-series solution for the fluid 
and metric variables in terms of t, and this convenient choice simplifies the 
equations. 

To obtain a correct power-series for the super-horizon initial condi- 
tions, we conduct a nonnal mode analysis. We begin with synchronous- 
gauge fluid and metric va riables JMa & Bertschinger 1995j ; iBucher et alj 
llOOCtlShaw & Lewisl201(]|) and equations of motion. We then define x = krj 
and new fluid variables, dividing out the relative factors of x. That is, we set 

(Al) 

(A2) 

(A3) 

(A4) 



S, = Si/x, 

ti = e,/x\ 

S'v = <^v/x, 
(7(3). 



Fl^'/.x", 



c-0) 



where Sj, Bj, Cy and Fj,'' are the usual fluid variables (e.g., see 
iMa & Bertschingejl 19951) . We then form a vector of the synchronous-gauge 
fluid and metric vaiiables, UZ = \6y,6y,6c,Sh,tyh,tvJc,'^v,f^,®,Tlm), 
where for this analysis we work at times early enough that the tightly 
coupled-photon baiyon fluid has a single velocity. In terms of the 
synchronous-gauge metric variable /;,„, we have = h[^, where ' denotes a 
derivative with respect to x. The full system of ODEs may then be written 



d\nx 



..A,/") Vk, 



(A5) 



where A^ are matrices containing coefficients of terms of different order 
in x; the matrices are obtained by Taylor-expanding the conformal Hubble 
parameter "K = aja = a)(2T + 1) / [t(t + 1)], all homogeneous densities 
(baryons, CDM, neutrinos, and photons), and pressures in powers of r <s: 1 . 
The space of solutions is spanned (to lowest order) by the eigenvectors J/'i 
(with eigenvalue A) of AqI 



U^T) = J]c,x'U^ 



(A6) 



Here c,\ are coefficients setting the contribution of each noiTnal mode to the 
solution, and can be chosen so that fluid variables match initial conditions. 
The physical growing normal modes are [3 
• Adiabatic mode (AD): A= I, 

jjW Mill 1 (23 + 4R,) 



3' 3' 4' 4' 36' 36(15 +4R,)' 

2 4 



0, : 



,1,0 



3(15 + 4R,)' 21(15 + 4R^) 
Baryon isocurvature (BI) mode: A = -I, 

Ui-l). BI ^ ^Q^ Q Q J^ ^ Q^ Q^ Q^ Q^ Q^ Q^T ^ 

CDM isocurvature (CI) mode: A = -I, 

{/!""• ° = (0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0)'' , 

k 

Neutrino density isocuivature (NDI) modes: A = -\, 
J/1""- ^°' = [-Ry/Ry, 1, 0, 0, -RyliARy), 1/4, 0, 0, 0, 0)^ , (AlO) 



(A7) 
(A8) 
(A9) 



Neutrino velocity isocurvature (NVI) mode, A = -2, 
[/!,"-' = (0, 0, 0, 0, ~Ry/Ry, 1, 0, 0, 0, 0, o)^ 



(All) 



Around each normal mode, we can extend to a solution liiM that 
includes higher order corrections: 



tli'kr) = ui'V + C/W x^+' + ... + [/!'' 



*,(l) 



kXi) 



+ ...., 



(A12) 



Modes with A <0 may still be 'growing' modes, ii ths physical variables 
(5, 6, o", ... oc jT^ for 7 > 0. 
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where the label Ui denotes the i'"" -order correction. We derive the cor- 

rections to the lowe st-order solution b y applying Eq. jASt to the ansatz, 
Eq UAIS . obtaining )Doran et al.l2003l) : 



(-1) 



[u+i)j-Ao][/<;;;^=A.y 

[(^ + 3)J-A„1{7!,-" =A,J/!,'" +A,f/';" +A.U 

I —"I t,(3) — ' t,(2) — - *,(1) — ^ ; 

1 -"J t,(4) — ' i.(3) -- A-.(2) 

+ A,[/':" +A4C/W. 






(A13a) 
(A13b) 
(A13c) 

(A 13d) 



Here I is the identity matrix in the space of all fluid+metric vaiiables. The 
solutions to this linear system can yield higher-order corrections to the time- 
evolution of the fluid variables for each normal mode. For density isocur- 
vature modes, the linear system Eq. l lA13bt becomes under-constrained. By 
directly applying 



l^r, 



"ir 

2 



~47TGa^Sp, 



(A14) 



Al Adiabatic mode (AD) 

The initial conditions for the well-known adiabatic mode are 
10 2 



(5. 

<5c = i5b 



(15 


+ 4Ry) 


3a,, 


(10 


+ 4Ry) 


byt// 


(15 


+ 4Ry) 




20 





a5 + 4Ry) 
3 

-6y = 30y.O 

1 



-2tfi 



e^ = e. 



(Ty 



Kkr^f 






(A16) 



3(15+4R,) 
It directly follows that the total initial entropy perturbation vanishes: 

S(0, k) = R^S, + RbSb - (3i4){RySy + RySy) 

= 6, + RbiSb - <5c) - (3/4)[(5y + Ry(dy - <5^)] ^ (A17) 

Thus this is an adiabatic (isentropic) initial condition, and initially only cur- 
vature perturbations are present. A curvature perturbations of amplitude 
f (0, k) causes a potential perturbation 1// = 2f (0, k)/(3ay). Therefore, adi- 
abatic modes enter the hoiizon with WKB amplitude A ^ -(3/2)t//(0,k) 
fau&Sugivama_1996:) . 



the Einstein constraint equation, however, we may close the system to ob- 
tain Ul and continue to higher orders using Eqs. l lA13ct - fA13dt . We 
th us reproduce the po wer series for the adiabatic and isocurvature modes 
of lBucheret^ i2000l) and lshaw & Lewis! ilOld) . 

As a final step, we must perform a gauge transformation to CN gauge, 
in order to have initial conditions for CosmoTherm. The gauge transform is 
given by 



(A15a) 

(A 15b) 

(A15c) 
(A15d) 

(A15e) 

(A15f) 
(A15g) 



^co„ 


= 5^ + 0-^, 




P, 


^o„ 


= ff; + ak^. 


^.pcon 


= 6P'; + aJ, 


0-™" 


= irl 


vrP t1 


i'm + dflm 



(p =rim -1-1 a. 



where superscripts 's' denote synchronous gauge variables, 'con' denote 
CN-gauge variables, 6Pi is the pressure perturbation of the f* species, 
and o", the anisotropic stress in the ('"' species. The homogeneous den- 
sity/pressure of the i'"" species are marked p, and P,, while dots denote 
derivatives with respect to conformal time. The results of the gauge trans- 
formation are listed below. 

CosmoTherm evolves temperature variables for the photons and neu- 



trinos, 0f, so that for example 6y = 40 



y.O. 



dy 



3^:01, and (Ty = 20 



r,2- 



The stiff ODE solver of CosmoTherm can readily compute their evolution 
correctly, even if the only the leading order terms for the initial conditions 
are included. That is, starting the calculation at very early times, well within 
the super-horizon regime, terms of order (kr^) in the initial conditions can 
in principle be neglected without significantly affecting the solutions. 



A2 Baryon and CDM isocurvature modes (BI/CI) 

For baryon isocurvature modes we have 
R^r (\5-4Ry) 



>P-- 



6y = Sy 

<5c = <5b 



2 (15 + 2R,) 

RbT (\5 + 4Ry) 

2 (15 + 2R,) 

(15+4R,) 



-^ 



-2i?br 



(15 + 2R,) 



3^ (15+4R,,) 
- - RuT 

2 (15 + 2R,) 



:3<A 



Oy = Oy 



^^ = -y«^^i5 



k^r] 



- 2Ry) 



RbT (15 
6 



4Ry) 



(15 + 2Ry) 

2RbT{krif 4i/r 



k-n= i;ipk r] 



(kr,r 



'I' ,2 
-kn 



:6yn. 



(A18) 



3(15 + 2R,) 45cv ' 45 

The initial conditions for CDM isocurvature modes can be obtained from 
these expressions by replacing R\y — > R^ and setting 

- 4Ry) 



■.s,-\ 



2 (15- 



2Ry) 



(A19) 



The baryon and CDM isocurvature modes both have vanishing initial po- 
tential perturbations (isocurvature condition) but total entropy perturbation 
S(i),k) ^ Ri + 0(t) for i e {b,c|. This means that a baiyon/CDM den- 
sity perturbation with amplitude di(0,k) leads to an entropy perturbation 
Si(0,k) = (Q.[/Q.iji)d[(0,k) initially. We also co mpute the gauge-inva riant 
Bardeen curvature variable, (, using the relation tShaw & Lewisll20ld) 



i-- 



. + 2 



> + '^/'H) 



(A20) 



3 ( 1 + w) 

where w is the total cosmic equation of state (equal to 1/3 during radi- 
ation domination). Using Eqs. JA18K it is easy to see that f = for 
the BI/CI modes when r = r; = 0, as should be the case for isocurva- 
ture perturbations. These initial conditions only excite the sin(/:rs) term 
of the photon monopole transfer function. We therefore have A = and 
B ^ -V6/4(ni/nm)(A:cci /«:)(! " 4 Ry/\5)6i(0,k), where again Ry accounts 
for the effect of anisotropic stress l lHu & Sugivamall 19960 . 
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A3 Neutrino density isocurvature mode (NDI) 

Although these modes are called isocurvature modes, they do not have van- 
ishing initial potential perturbations: 

2R,, _ IRy 

(15+4R,,) ~"T5fl^ 
Rv 1 , 



^ = ~-, 



(15+4i?,) 2' 
Ry{n + 8R,,) __Rv 

R^(15 + ARy) ~ Ry 

(15 + 8i?,) 



Sc= S\y- 



(15 + ARy) 

3Ry 



l+< 



0y = 0b 



(l5+4Ry) 

\9Ryk^r] 

~ ' ARy(\5 + ARy) 

l5k^T] 15 
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^^•^f;: 



2i?, 



■i/r^:-;; 



4(15 + 4R^) 



ipk^Tj 



Ry 



■l?n= -il/k^ri 
(15+4R,,) ' 2^ ' 



{krif 



- 2(15. 4«.) ■-iy^'''^" = Wy^"^^" = T-5'"^- (^2>) 

An initial neutrino density perturbation of (5,,(0, k) leads to entropy pertur- 
bation S (0, A:) ^A(l>(0,k) ^ R,,(5y(0,j(:)/(15 + 4R;,) and potential perturbation 
t//{0,k) =! -2RySy{0,k)/{l5ay). Correspondingly, these modes excite both 
sin(irs) and cos(<:rs) parts of the photon monopole transfer function. As for 
the BI/CI modes, ( = initially for the NDI mode. 



A4 Neutrino velocity isocurvature mode (NVI) 

The second isocurvature mode for neutrinos is sourced by non-vanishing 
velocity perturbations relative to the photons: 

4Ry 



If,-- 



(5 + 4-Ry)krj 

4Ry 

(5 + ARy)kn 
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~ Ry{5 + ARy) ARy ^ '^ 



5k 
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ARyk 
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Akr] 



k + if/k-^T] 
= 4,k^r, 



3R, 



(kilY 



(A22) 



3(5 + ARy) 

These modes are very hard to exc ite and also show a coordinate divergence 
for small krj JBucher et al.ll2000l) . Numerically it is straightforward to in- 
tegrate the corresponding perturbation equations with CosmoTherm. Also, 
all physical quantities are gauge-independent and we find these to converge 
very well. As for the NDI mode, f = initially for the NVI mode. 
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